Optimal fluctuation approach to a directed polymer in a random medium 
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A modification of the optimal fluctuation approach is applied to study the tails of the free energy 
distribution function Pl(F) for an elastic string in quenched disorder both in the regions of the 
universal behavior of Pl (F) and in the regions of large fluctuations, where the behavior of Pl (F) is 
non-universal. The difference between the two regimes is shown to consist in whether it is necessary 
or not to take into account the renormalization of parameters by the fluctuations of disorder in the 
vicinity of the optimal fluctuation. 
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A large variety of physical systems can be described in 
terms of an elastic string interacting with a quenched ran- 
dom potential. The role of such a string can be played by 
a domain wall in a two-dimensional magnet, a vortex line 
in a superconductor, a dislocation in a crystal and so on, 
however most often this class of systems is discussed un- 
der the name of a directed polymer in a random medium. 
The unfading interest to this problem is additionally sup- 
ported by its resemblance to more complex systems with 
quenched disorder (e.g., spin glasses), as well as by its 
relation to the dynamics of a randomly stirred fluid and 
to stochastic growth (see Refs. [l] and@ for reviews). 

A number of properties of a random directed poly- 
mer can be easily deduced as soon as one knows Pl(F), 
the free energy distribution function for large polymer 
length L. In 1987 Kardar— proposed a method for cal- 
culating the moments of Pl{F) in the one-dimensional 
(ID) problem with (5-correlated potential. However, later 
it was understood^ that the approach of Ref. 3 allows one 
to find^ only the tail of Pl(F) at large negative F (the 
left tail), so the conclusions on the width of this distri- 
bution function have to rely on the assumption that at 
large L it acquires a universal form, 



P L (F)=P*(F/F*)/F*, 



(1) 



incorporating the dependence on all parameters through 
a single characteristic free-energy scale F*(L) oc L u , 
which therefore can be extracted from the form of the 
tail. Here and below it is assumed that F is counted off 
from its average. 

In the present work we report the results of the first 
systematic investigation of Pl (F) in the directed polymer 
problem which makes a distinction between the interval 
of the universal behavior and the regions of large fluc- 
tuations, where Pl(F) can be expected to deviate from 
Eq. (P). We demonstrate that the far- left tail of P L (F) 
is determined by the optimal fluctuation and apply an 
analogous approach for the analysis of the far-right tail 6 . 
Inside the universal region, we rely on the generalization 
of this method accounting for the renormalization of dif- 
ferent parameters by fluctuations. 

This allows us to show that in the ID case the exponent 
r)+ in the dependence In Pl(F) oc — \F\ n + is changed from 
3 in the universal part of the right tail to 5/2 outside of it 



(at the largest F). On the other hand, we find that in the 
left tail the behavior of Pl (F) is qualitatively the same 
(with r)_ = 3/2) in both regions and that the analysis of 
Ref. reproduces a correct estimate for F * (L) only as a 
result of this coincidence, because in reality it describes 
the non-universal part of the left tail. We also analyze 
how the tails of Pl (F) are modified when the distribution 
of disorder is characterized by a finite correlation length. 

The model. - We consider an elastic string confined to 
a plane and interacting with a random potential V(t, x) 
with zero mean and Gaussian statistics described by 



V{t, x)V(t', x 1 ) = 5(t - t')U(x - x') . 

An overbar denotes the average with respect to disorder 
and t is the coordinate along the string. Such a string 
can be described by the Hamiltonian, 





'dx(t')' 




dt' 



+ V[t',x(t')} 



(2) 



The partition function of a string which starts at t = 
and ends at the point (f, x) is denoted z(t, x) and its free 
energy f(t,x) = —Tlnz(t,x), T being the temperature. 

Derivation of f(t, x) with respect to t shows^ that the 
evolution of f(t, x) with the increase in t is governed by 
the Kardar-Parisi-Zhang (KPZ) equation 9 (with the in- 
verted sign of the non-linear term), 



ft + Q-/2J)%-vf xx = V(t,x) 



(3) 



where t plays the role of time, subscripts t and x de- 
note partial derivatives and v = T/2J. On the other 
hand, derivation of Eq. with respect to x leads^ to 
the Burgers equation* with random force, 



tit + uu x - vu xx = V x (t, x)/J , 



(4) 



in which u(t, x) = f x (t, x)j J plays the role of velocity. 

The replica approach to the directed polymer problem 
is based on calculating the moments, Z n = Z n , of the dis- 
tribution of the partition function, Z = z(L, 0), and does 
not rely on the analytical continuation of n to 0. Kardar— 
has noticed that for any integer n > 1 (and large enough 
polymer length L) Z n with an exponential accuracy can 



2 



be approximated by exp[— Eq(ti)L/T\, where Eq(ji) is the 
ground state energy of n bosons whose mass is equal to 
J/T 2 and interaction to —U(x)/T. On the other hand, 
since Z n = exp(—nF/T), Z n can also be expressed in 
terms of Pl{F), the distribution function of the free en- 
ergyF = /(L,0) = -TlnZ: 



dFP L (F) exp(-nF/T) . 



(5) 



In the case of (5-correlated random potential, 
U(x) = Uo5(x), Eo(n) can be found exactly^. This gives 
Z n oc exp[(JU 2 /2AT 5 )n 3 L], where the linear in n term 
in Eo(n) has been omitted, because it can be elimi- 
nated by a constant shift of the potential V(t, x) in 
Eq. ([2]). Comparison of the two expressions for Z n al- 
lows one to conclude^ that for large negative fluctuations 
P L (F) oc exp[-S(F)], where 



3/2 



JU§L 

J^2 



1/3 



(0) 



For such Pl(F) the integral in Eq. ([5]) at positive n is 
dominated by the saddle point at F « —F c (L)n 2 , where 

F C (L) = F 3 (L)/T 2 = JU 2 L/8T 4 . (7) 

Since the explicit expression for Z n is known only 
at integer n > 1 and cannot be analytically continued 
to < n < the region of the applicability of Eq. 
Q is restricted to —F^>F C (L). With increase in L 
this region is shifted to larger and larger values of 
the ratio \F\/F 9f 3> F c /F* oc L 2 / 3 , which shows that the 
replica approach cannot provide any information on 
P*[F/F*(L)], *- e - on the universal form of Pl(F) at 
L — * oo. 

The optimal fluctuation approach. - If V(t, x) would 
be independent of t, for large enough L the free energy 
of a polymer in a given disorder realization would be 
given by F « EqL, where Eq is the ground state energy 
of the single-particle quantum-mechanical Hamiltonian, 
H = -(T 2 /2J)d 2 + V(x). When V is fluctuating with t, 
one can nonetheless consider a question about the form 
of the most probable "stationary" (i.e., uniform along t) 
fluctuation of V which leads to the given negative value 
of Eq. According to the analysis of Ref. [ID, this form is 
given by 



V(x) = -A 



dx'U(x-x')^> 2 (x'), 



(8) 



where A is an arbitrary constant and ^(x) is the lo- 
calized solution of the non-linear Schrodinger equation, 
Eq^/ = , in which V(x) has to be replaced by the 
right-hand side of Eq. . 

In the case of (5-correlated random potential, 
U(x) = Uo5(x), such a solution (the soliton) has a formi 2 - 
^a(x) = T/[(AJU y/ 2 A cosh(>/A)], where the param- 
eter A(F) = T[L/2J(-F)} 1 / 2 can be called soliton 
width. Substitution of ^a(^) into Eq. ([H]) then gives 

V(x) = -{T 2 /JA 2 ) cosir 2 (x/A) . (9) 



In such a potential there exists only one level with a 
negative energy 1 - , from where the condition for the ap- 
plicability of the relation F « E L is —F^>T. 

On the other hand, substitution of Eq. ^ into the 
functional, 



5' 



1 

2t/o 



dt 



+ 00 



dxV 2 (t,x) , 



(10) 



determining the probability of V(t,x) for U(x) = Uo5(x) 
leads to exactly the same expression for S(F) as follows 
from the replica analysis, see Eq. ([6]). This proves that for 
5-correlated random potential the far- left tail of Pl (F) is 
entirely determined by the optimal fluctuation of V(t, x). 
However, in order to find the range of the applicability 
of such an approach one has to consider the influence of 
the fluctuations of disorder in the vicinity of the optimal 
fluctuation. This can be done most transparently in terms 
of the Burgers equation representation. 

The Burgers equation approach. - It follows from the 
definition of u(t,x) = f x (t,x)/J that the stationary so- 
lution of the Burgers equation (j4|) with the pumping po- 
tential V(x) of the form ((9J is given by 

T d 

ua(x) = — -— ln^ A (a;) = u A tanhfz/A) , (11) 
J ax 

where ua = T/JA. Thus, in terms of the velocity u(t, x) 
the stationary soliton looks like an inverted standing 
shock wave^ Such solitons have been discussed (in var- 
ious contexts) in a number of works by Fogedby^ 

The description in terms of the Burgers equation allows 
one to understand rather easily how the shape of the 
optimal fluctuation evolves with time when one takes into 
account the initial condition. Below we adopt the free 
initial condition for the string, z(0, x) = 1, which in terms 
of interface dynamics corresponds to the standard initial 
condition 2 , for non-stationary growth, u(0,x) = 0. 

If one assumes that for t > the potential V(x) is of 
the form |9]), then the solution of (j4|) has to be close to 
(fTTj) for not too large x, but must be close to zero at large 
enough distance from the region around x = where the 
pumping is localized. It is clear that this can be realized 
via the presence of two travelling shock waves, which docs 
not require any additional pumping and therefore does 
not change S(F). The velocities of these waves are deter- 
mined by their amplitude, v = ±ua/2, so for t A/ua 
they will be located at \x\ ~ UAt/2 ^> A, i.e., relatively 
far from the soliton core. Accordingly, in this regime 
F = J J ( ^ oo dxu(x) has to be close to —(J/2)u A t, as it 
has been assumed above. It is possible to check that the 
imposition of the fixed initial condition [x(t — 0) = 0] 
also does not lead to the change of S(F). 

In terms of the Burgers equation parameters (the 
viscosity v = T/2J and the pumping force intensity 
D = Uq/2J 2 ), Eq. ([6|) can be rewritten as 



^ - Is 



2{-Ff 
J 3 L 



1/2 



(12) 
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From the nature of the optimal fluctuation approach it is 
clear that this expression can be expected to be valid 
only when the soliton is so narrow, A <c xq (where 
xq ~ v 3 / D ~ T 3 / JUq), that one can neglect the renor- 
malization of parameters by the nonlinearity. Note that 
the constraint A <C xq is equivalent to the condition 
— F 3> F C (L) defining the region of the applicability of 
the replica approach results. 

At — F <C F C (L) (i.e., A 3> x ) one has to take into 
account the influence of the fluctuations in the vicinity 
of optimal fluctuation. Since optimal fluctuation is quasi- 
stationary (see below), this (up to a numerical factor) can 
be done by replacing all parameters in Eq. (fi"2]) by their 
effective values at the scales of the order of A (and zero 
frequency). However, it is well known that the amplitude 
of the non-linear term in Eq. ^ (and, therefore, coef- 
ficient J) cannot be renormalized-f due to the Galilean 
invariance. On the other hand, in the ID case the ratio 
H = D jv has to remain unchanged 9 - as a consequence of 
the fluctuation-dissipation theorem obeyed by Eq. (g]). 8 
The coefficient /i describes the spectral density of the 
equilibrium equal-time fluctuations of u, which is known 
to be exactly the same as in the absence of the non- 
linearity4 

Since Eq. (TT2")) includes only invariant coefficients 
H = D jv and J, it can be expected to remain appli- 
cable not only for — F 3> F C (L), but also (with some 
numerical correction) in the region of the universal be- 
havior, —F < F C (L). From the other side the region of 
the applicability of this estimate is restricted by the con- 
straint L (A 3 //!) 1 / 2 , which is required—, for the at- 
tainment of equilibrium at the length-scales of the or- 
der of A and is satisfied as soon as — F ^> F*(L) [i.e., 
S(F) 3> 1]. The fulfillment of this condition is neces- 
sary for the applicability of the above argument (based 
on known properties'^ of an equilibrated system) in our 
non-stationary situation. 

The right tail. - The very special shape [in terms 
of fit, x)] of the optimal fluctuation corresponding to 
the left tail originates from the possibility to have a 
growing fluctuation inside of which almost everywhere 
ft + (1/2 J) f 2 as 0, and therefore only a small part of the 
volume of this fluctuation contributes to the functional 
(fTUJ) [where V(t, x) should be replaced by the left-hand 
side of Eq. ([3])]. It is clear that in the case of the right 
tail such a situation is impossible because ft has to have 
the same sign as /J. As a consequence, the optimal fluc- 
tuation corresponding to the right tail must have a more 
simple shape which can be characterized by a single rele- 
vant length-scale, A. This length-scale can be estimated 
from ft ~ (1/2 J) fl, which gives A(F) ~ (FL/J) 1 / 2 and 
a temperature-independent estimate for S(F), 



S(F) 



i?5/2 



C/ J 1/2 i 1/2 



(13) 



Note that Eq. (TT3)) contains the unrenormalized coeffi- 
cient Uq, whereas the estimate for A(F) shows that it 
grows with increase in F. This suggests that even for 



large F it may be necessary to correct Eq. (TT3l) by taking 
into account the renormalization of Uq. 

It follows from the analysis of higher-order diagrams 18 
that in a stationary ID system the scale dependence 
of v and D = Uq/2J 2 can be found from the self- 
consistent one-loop theor y 19 ' 20 (mode coupling approxi- 
mation), which for A 3> xq give o 17 ' 20 v{A) ~ (^A) 1 / 2 , 
D(A) ~ (/x 3 A) 1 / 2 . However, in the case of the right tail 
the replacement of D by D(A) (where A is the optimal 
fluctuation width) overestimates the renormalization of 
D. It is so because the only source for the renormaliza- 
tion of v and D is the non-linear term in Eq. and 
therefore this renormalization can be effective only up to 
a length-scale R ~ fiJL/F at which the typical veloc- 
ity of equilibrium fluctuations, ur ~ (/i/i?) 1 / 2 , becomes 
comparable with the velocity up ~ F/JA(F) character- 
izing the optimal fluctuation being considered. For the 
left tail R ~ A, whereas for the right tail R <C A, so in 
the latter case the effective value of D has to be estimated 
as D(R) ~ ^(JL/F) 1 / 2 < D(A), from where 



S(F) 



DiR)^/ 2 ^/ 2 



(14) 



The conditions for the applicability of Eq. (| 14[) describ- 
ing the universal part of the right tail are S(F) ^> 1 
and R > x , which means F*(L) < F < F C (L). 
Note that here both F*(L) and F C (L) are the same as 
in the left tail, i.e. are given by Eqs. ([6]) and ((7]) re- 
spectively. Another similarity to the left tail is that the 
fulfillment of the condition F ^> F*(L) ensures that 
L > A 2 /^(A) ~ (A 3 /^) 1/2 , and therefore time L is suf- 
ficient for attaining the equilibrium at the length-scales 
of the order of AJZ 

At F ^> F C (L) one gets R <C xq, which means that the 
renormalization of v and D becomes completely ineffec- 
tive. Therefore, the far-right tail is described by Eq. (jT3j) 
with unrenormalized parameters. 

A finite correlation length. - If correlations of ran- 
dom potential are characterized by a finite correlation 
length £ = Uq/U(0) <C xq [where now Uq = J dxU(x)], 
the expression ^ remains applicable only for A 3> £, 
that is for — F <C F^ ~ (xq/^) 2 F c . In the other regime 
(at — F ^> F%) the characteristic size of Vt'(x) becomes 
much smaller then £, which leads 2 ^ to V(x) oc —U(x), 
E cx -V(0) and S « F 2 /2U(0)L. This means that the 
finiteness of £ makes the most distant part of the non- 
equilibrium left tail Gaussian 22 and independent of T. 

With decrease of xq [i.e., of temperature) the region 
with Gaussian behavior becomes larger and larger, and 
for xq < f [i.e., T < T ~ (J[/ £) 1/3 ] the whole 
non-universal tail must be Gaussian. In this regime 
the parameters of the universal part of the left tail 
also are modified. For x -C £ the spectral density of 
the equilibrium fluctuations of u(x) at scales exceed- 
ing £ can be estimated as where ~ (Dt^/^ 3 ) 1 / 2 
is the characteristic velocity which is created by ran- 
dom force with characteristic length-scale £ during the 
time ~ £/ M £ required for the breaking of such a 
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fluctuation. This gives one an estimate fi ~ (D 2 /C) 1 ^ 3 , 
substitution of which into Eq. (fT2"|) instead of the ra- 
tio D/v leads to S(F) ~ [T/T*] 3 / 2 , where the free 
energy scale F*(L) ~ (x n /£) 2 / 9 F*(L) ~ ( J[/ 4 /e 2 ) 1/9 i 1/3 
(replacing T*) does not depend on T. Directly in terms of 
the polymer problem an analogous estimate for T* (L) can 
be obtained from scaling arguments complemented by the 
assumption that at low enough temperatures T*(T) has 
to be temperature independent r 3 - and follows also from 
the replica-symmetry breaking analysis of Ref. [2J. 

It can be shown that F*(L) replaces F*(L) also in 
Eq. (fT4"]) describing the universal part of the right tail. 
On the other hand, the fmiteness of £ does not change 
the form of the non-universal part of the right tail de- 
scribed by Eq. (fl~3|) . All this allows one to conclude 
that at low temperatures Pl{F) becomes temperature 
independent, whereas the crossover between the regions 
of universal and non-universal behavior takes place at 
\F\ ~F C ~ (U 2 /Je) 1/3 L. 

Conclusion. - In the present work we have studied the 
form of the free energy distribution function Pl(F) in 
the ID directed polymer problem both for ^-correlated 
random potential and for a finite correlation length £. In 
terms of the KPZ problem the same distribution function 
describes the distribution of heights in the process of non- 
stationary growth starting from flat interface. 

Our analysis has shown that for \F\ <C F C (L) oc L both 



tails of Pl (F) are characterized by the same free energy 
scale, F*(L), which confirms the universality hypothesis 
((U). For T » T ~ (J[/o£) 1/3 tnis sca l e is temperature 
dependent [see Eq. ©] , whereas at T <C To it saturates at 
a finite value. On the other hand, the exponents n- =3/2 
and r] + = 3 in the dependence In Pl(±F) oc~(F/F r ) v± 
are insensitive to the relation between T and Tq. The 
same values of ry_ and n + have been found by Prahofer 
and Spohii 2 ^ for the polynuclear growth (PNG) model, 
which in terms of a directed polymer problem, Eq. (|2j). 
corresponds to the Poisson distribution of idential point- 
like impurities and a rather perculiar limit of J = and 
T = . 25 i 26 This agreement confirms that the two models 
indeed belong to the same universality class. 

We also have investigated the non-universal tails of 
Pl(F) (at \F\ 3> F c ). In particular, we have demon- 
strated that in the far-right tail rj + = 5/2, which is in 
agreement with the numerical results of Kim et al.^ 
who found n + = 2.4 ± 0.2. On the other hand, the form 
of the far-left tail depends on £. For £ = it is the same 
(?7_ = 3/2) as in the universal region, whereas for £ > 
the left tail becomes Gaussian for sufficiently large \F\. 

The authors are grateful to V.B. Gcshkcnbcin, A.I. 
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sions. The work of I.V.K. was supported by RFBR under 
grant 05-02-17305. 



1 M. Kardar, Lectures on Directed Paths in Random Media 
(Les Houches, 1994), |cond- mat/941 1022] 

2 T. Halpin-Healy and Y.-C. Zhang, Phys. Rep. 254, 215 
(1995). 

3 M. Kardar, Nucl. Phys. B 290 [FS 20], 582 (1987). 

4 E. Medina and M. Kardar, J. Stat. Phys. 71, 967 (1993). 

5 Y.-C. Zhang, Europhys. Lett. 9, 113 (1989). 

6 An analogous approach has been used^ to study the dis- 
tribution of of velocity and its derivatives in the Burgers 
turbulence problem, which however requires to consider 
optimal fluctuations with completely different structures 
then here. 

7 E. Balkovsky, G. Falkovich, I. Kolokolov and V. Lebedev, 
Phys. Rev. Lett. 78, 1452 (1997); A.I. Chernykh and M.G. 
Stepanov, Phys. Rev. E 64, 026306 (2001), see also M.V. 
Feigel'man, Zh. Eksp. Teor. Fiz. 79, 1095 (1980) [Sov. 
Phys. - JETP 52, 555 (1980)]. 

8 D.A. Huse, C.L. Henley and D.S. Fisher, Phys. Rev. Lett. 

55, 2924 (1985). 

9 M. Kardar, G. Parisi and Y.-C. Zhang, Phys. Rev. Lett. 

56, 889 (1986). 

J.M. Burgers, The Nonlinear Diffusion Equation (Reidel, 
Boston, 1974). 

1 H.B. Thacker, Rev. Mod. Phys. 53, 253 (1981). 

2 B.I. Halperin and M. Lax, Phys. Rev. 148, 722 (1966). 

3 L.D. Landau and E.M. Lifshitz, Quantum Mechanics 
(Nauka, Moscow, 2002; Pergamon, New York, 1981), Sec. 
23, Problem 5. 

4 This means that u(x) = —ua(x) describes the genuine 
standing shock wave which for any A is the solution of Eq. 



Q with the vanishing right-hand side. 
H.C. Fogedby, Phys. Rev. E 68, 026132 (2003) and refer- 
ences within. 

E. Medina, T. Hwa, M. Kardar and Y.-C. Zhang, Phys. 
Rev. A 39, 3053 (1989). 

T. Nattermann and L.-H. Tang, Phys. Rev. A 45, 7156 
(1992). 

H.K. Janssen and B. Schmittmann, Z. Phys. B 63, 517 
(1986); V.S. L'vov, V.V. Lebedev, M. Paton and I. Pro- 
caccia, Nonlinearity 6, 25 (1993). 

H. van Beijeren, R. Kutner and H. Spohn, Phys. Rev. Lett. 
54, 2026 (1985). 

T. Hwa and E. Frey, Phys. Rev. A 44, R7873 (1991). 
B.I. Shklovskii and A.L. Efros, Electron Properties of 
Doped Semiconductors (Nauka, Moskow, 1979; Springer- 
Verlag, Berlin, 1982). 

This can be also shown in terms of the replica approach. 
The tail becomes Gaussian when the size of the ground- 
state wave-function of n bosons becomes much smaller 
than their interaction range £. 

T. Nattermann and W. Renz, Phys. Rev. B 38, 5184 
(1988). 

S.E. Korshunov and V.S. Dotsenko, J. Phys. A 31, 2591 
(1998). 

M. Prahofer and H. Spohn, Phys. Rev. Lett. 84, 4882 
(2000). 

K. Johansson, Comm. Math. Phys. 209, 437 (2000). 
J.M. Kim, M.A. Moore and A.J. Bray, Phys. Rev. A 44, 
2345 (1991). 



